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Abstract 

In this paper, we present an analysis of the outage probability for fixed-gain amplify-and-forward 
(AF) multihop relay links operating in the high SNR regime. Our analysis exploits properties of Mellin 
transforms to derive an asymptotic approximation that is accurate even when the per-hop channel gains 
adhere to completely different fading models. The main result contained in the paper is a general 
^vq ■ expression for the outage probability, which is a functional of the Mellin transforms of the per-hop 
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channel gains. Furthermore, we explicitly calculate the asymptotic outage probability for four different 
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systems, whereby in each system the per-hop channels adhere to either a Nakagami-m, Weibull, Rician, 
or Hoyt fading profile, but where the distributional parameters may differ from hop to hop. This analysis 
leads to our second main result, which is a semi-general closed-form formula for the outage probability 
of general fixed-gain AF multihop systems. We exploit this formula to analyze an example scenario 
for a four-hop system where the per-hop channels follow the four aforementioned fading models, i.e., 
the first channel is Nakagami-m fading, the second is Weibull fading, and so on. Finally, we provide 
simulation results to corroborate our analysis. 
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I. Introduction 

Amplify-and-forward (AF) relay systems have received a lot of attention recently due to their 
ability to improve coverage and, thus, capacity in a geographical sense. To date, two main AF 
protocols have been focused on in the literature: variable-gain (a.k.a. channel state information 
(CSI) assisted) AF relaying and fixed-gain AF relaying (see, e.g., [1J and the references therein). 
While the former method yields good performance when CSI is available at the relay nodes, 
the latter technique is more suitable in simple systems where such information is not available, 
although the performance of the system often suffers. In particular, fixed-gain AF relaying may 
be a good choice in low-complexity systems, such as emerging energy and utility management 
applications (e.g., "smart grid" and water metering communication networks) as well as industrial 
wireless sensor networks [|2]. Consequently, we focus on the fixed-gain protocol in this paper. 

As with many wireless communication systems operating in fading environments, the end-to- 
end outage probability is an important metric that can be used to characterize the performance 
of a fixed-gain AF relaying systemj- Several results on this topic have been published. In [OQ, 
the authors derived a bound on the end-to-end SNR of a fixed-gain AF link, which was used to 
study the outage probability when each hop fades according to a Nakagami-m distribution. This 
analysis was adapted and extended in [|3l for cases where the per-hop fading distributions are 
Nakagami-n (Rice) and Nakagami-g (Hoyt), and the outage probability was studied using Pade 
approximants. In flU, the authors derived tight closed-form bounds on the outage probability at 
asymptotically high SNR for the case where the underlying channel power probability density 
function (PDF) is nonzero at the origin, a condition that is valid for Rayleigh, Rician, and 
Hoyt fading, but excludes Nakagami-m (for m > 1) and Weibull fading. Other more recent 
performance analysis results for multihop relay networks focus on variable-gain AF relaying, 
particularly in Nakagami-m fading channels J6j, 0. 

There are two main drawbacks to results currently available on fixed-gain AF relaying, which 

'Symbol error probability is also an important metric that has been investigated for fixed-gain AF relaying (see, e.g., I3l-l5l 
and the references therein). 
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we aim to improve upon in this paper: 

1) Many of the results found in the literature are given as lower bounds on the outage 
probability (see, e.g., (D, 0, flU). To date, exact asymptotic results (not bounds) for 
general multihop systems have not been reported. Indeed, exact results for dual-hop systems 
operating in a class of Nakagami-m channels have only recently been published (8), ||5). 

2) Most analysis of multihop systems to date considers homogeneous fading, i.e., each hop 
fades according to the same distribution. Results for inhomogeneous systems appear to 
be limited to dual-hop links (see, e.g., ifTOl . ifTTlO . although this scenario is likely to be 
encountered frequently in practical multihop systems. 

In this paper, we present a general framework for analyzing the outage probability of fixed-gain 
AF multihop relay systems operating in the high SNR regime. In contrast to outage performance 
analysis that currently exists in the literature, our approach does not rely on bounding the 
outage probability, but rather exploits properties of Mellin transforms to derive an asymptotic 
approximation - which is in the form of a functional of the Mellin transforms of the per-hop 
channel gain PDFs - that is accurate even when the per-hop channel gains adhere to completely 
different fading models. Furthermore, the nature of the functional allows us to apply the residue 
theorem from complex analysis to derive asymptotic approximations for the outage probability 
for specific fading models - including Nakagami-m, Weibull, Rician, and Hoyt fading - in terms 
of elementary functions, which can be calculated easily in practice. Our contribution culminates 
in the introduction of a semi-general asymptotic formula for the outage probability of fixed-gain 
AF multihop systems. 

The rest of the paper is organized as follows. In Section II, we define the AF system model. In 
Section III, the general framework for the asymptotic outage probability is detailed. This analysis 
draws heavily on properties of Mellin transforms (outlined in Appendix A for convenience). In 
Section IV, the analytical framework is applied to scenarios whereby the per-hop channels adhere 
to either a Nakagami-m, Weibull, Rician, or Hoyt fading profile, or a combination thereof. This 
section concludes with the introduction of a semi-general asymptotic formula for the outage 
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probability and a brief discussion on the convergence of the asymptotics. Simulation results that 
corroborate our analysis are presented in Section V, and conclusions are drawn in Section VI. 



II. System Model 

Consider a multihop network with a source node, a destination node, and N — 1 relay nodes in 
between where N > 2 (see Fig. [TJ). Communication can only be achieved in a half-duplex manner 
between adjacent nodes. A data symbol d is conveyed to the first relay node. For simplicity and 
without loss of generality, we let E [|<i| 2 ] = 1. This symbol is affected by flat fading in the 
transmission medium and additive Gaussian noise at the receiver (i.e., the relay node). 

The received signal is then amplified by a fixed gain A\, then conveyed to the next relay 
node and so on until the destination is reached. Denote the channel coefficient modelling the 
channel between the (n — l)th relay (or the source) and the nth relay (or the destination) by 
h n . Also, denote the additive noise at the nth relay node (or the destination) by v n , which is 
zero-mean complex Gaussian distributed with variance N >n /2 per dimension. Now, we can write 
the following input-output system equation Ifl2ll 

N \ N / N \ 

yi A n _ x h n \d + n A i- ih 3 ) Vn (i) 

\n=l / n=l \j'=n+l / 

where A = 1. Various amplification factors have been proposed in the literature, including flTJI 

A n = . 1 n = l,...,N-l. (2) 



'E [\h n f]+N 0>n 

In order to maintain generality, however, we do not explicitly define A n in what follows. 

Since we will eventually be interested in the asymptotics of the outage probability, it is 
beneficial to define a reference parameter 7 that gives a notion of average SNR across the N- 
hop link. In this case, we let N 0jn = p n j^ for n — 1, . . . , N where {p n } are strictly positive (and 
finite) scaling factors and p 1 = 1 for convenience. In order to illustrate the physical definition 
of 7, we note that when © is adopted for the amplification factors, the average SNR for the 
nth hop is given by j n = (E [\h n \ 2 ~\ / p n ) 7. We can now write the following expression for the 
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instantaneous end-to-end SNR: 



rin=l IM 



SNR = N ii^i^zU _^ (3) 

Z] n =l Pn rij=„+l \hj I 

III. Outage Probability 

For the ease of exposition, we define the random variable X ri = A^_ x \h n \ 2 in the following 
analysis. It follows from © that the outage probability of the fixed-gain AF multihop link can 
be written as 

P = P (SNR < lth ) 

(N N-l N N 

n=l n=l j=n+l 

= P\^---((X l -a l )X 2 - a 2 ) ■■■^X N < a N j (4) 
where we have defined a n (7) = pnlth/l for brevity. We further define the random variables 

Z n = W n X n+1 , n = 0,...,N-l (5) 

and 

W n = W n . X X n - O n > 0, 71 = 1, . . . , N - 1 (6) 

with Wo = 1. Note that W n is a conditional random variable in that it relates to the translation 
of Z n _i where it is given that Z n _ x > a n . Also, it is clear that W n and X n+X are statistically 
independent. Now we can apply this same conditioning on © recursively to obtain 

Po = P (Z <a 1 ) + P(Z Q >a 1 )p[(--- (W X X 2 -a 2 )---)X N < a N ) 
= P(Z <a 1 ) + P(Z > a x ) [P (Z x <a 2 ) + P (Z x >a 2 ) ■ ■ ■ 

(P (Z N - 2 < a N „ x ) + P (Z N - 2 > (jjv-i) P (Z N _! < a N ) ) • • • ) 

= i - n ^ m . (7) 

n=l 
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This is a satisfyingly simple exact expression for the outage probability, although the calculation 
of the CCDFs of {Z n } remain. It turns out we can construct an elegant lemma using Mellin 
transforms (see Appendix A) that aids this calculation. 

Lemma 1: The Mellin transform of fz n can be approximated by 

M [fz n ; s] » 1 E • • • E (-l) E ^ J 

n+1 

[/^-js-EL/fc] (8) 



r( S )n; =1 ^ 

X 



'■(-s'-E^i)ii; ,^ 

where L%, . . . ,L n > 1 are integers that define the order of the approximation. 

Proof: First, we note that M [f Zo ; s] = M [f Xl ; s}. Now, from © and property (l32i we 
have that M [fz n+1 ]s\ = M \fw n +i > s ] ^ [/x„ +2 ; s] . It is easy to see from the definition of the 
random variable W n that fw n+1 (w) = fz n (w + <y n +i) / Fz n (cn+i) for u> > 0. The proof follows 
from induction on n, where properties (|29l l and (|30T ) are applied in the inductive step. ■ 

We are now in the position to state our first main result in the form of the following proposition. 

Proposition 2: The outage probability of an iV-hop fixed-gain AF link is asymptotically given 

by 

Li-l ijv-i-1 N-1 1 j ^e n 



li=0 £jv_i=0 n=l x/ 



1 r* + "° /'WttV l' (Aft,-! -£) fr 

as 7 — )■ oo with k > Ajv-i, where X n = E&=iA and fx n denotes the PDF of the channel power 
for the nth hopP. 

Proof: Starting from Lemma [U we use (1311 to calculate M \Tz N ^\ s], then take the inverse 
transform using (|28l ) and substitute this into ©, resulting in a cancellation of the product 
n^Ti 1 -^Zn-i ( ff n)- Eq. © follows from a change of variables in the contour integral. ■ 



2 The notation ~ denotes asymptotic equivalence in the relevant variable and limit, 7 — > 00 in this case. 
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Although this appears to be a rather complicated expression for P Q , we note that it can be 
computed as long as we know the Mellin transforms of the individual PDFs f Xn , which can 
generally be calculated easily for many cases of interest. Thus, assuming we can compute these 
transforms, we can at least evaluate the outage probability bound numerically to a high degree 
of accuracy provided the transforms decay quickly as |Im(s)| — > oo. 

The Mellin transforms for the main fading distributions of interest are given in Table [jf . The 
key thing to notice from this table that makes the ensuing analysis uniform and tractable is that 
each PDF decays exponentially. In the transform domain, this results in the property that the 
transform of each fading distribution has poles at various points along the real axis. We will 
exploit this property in the next section to derive simple asymptotic expressions for P Q under 
various fading conditions. 



IV. Analysis of Common Fading Distributions 

We now apply the general analysis detailed above to a number of common fading distributions. 
For much of this section, we assume that all hops adhere to the same class of distribution 
{e.g., Nakagami-m) but where the channels for individual hops may vary in their distributional 
parameters {e.g., shape factors). We label this condition "homogeneity" in this context, and we 
later relax this restriction in order to analyze general multihop links as well as to show the power 
and versatility of adopting the proposed analytical framework. The distributions considered here 
(also listed in Table U) are derived from Nakagami-m, Weibull, Rician, and Hoyt fading. The 
analysis can be extended to other fading distributions using © and the techniques outlined 
below. The section concludes with our second major result, which is in the form of a semi- 
general closed-form formula for the outage probability for large 7. 

3 The transforms for the Rician and Hoyt distributions can be calculated by consulting standard tables of integrals and transforms 
{e.g., d, El). 
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A. Nakagami-m Fading 

If the per hop channels adhere to a Nakagami-m fading profile, the random variable X n has 
a gamma density function with scale parameter 8 n and shape parameter m n . The corresponding 
PDF and Mellin transform are given in Table |U By substituting this transform into © and 
applying the Mellin-Barnes integral definition of the Meijer G-function |fT5l along with the 
functional relations ITT31 9.31.2] and [[T51 9.31.5], we can write the outage probability as 



Ajv_i c n N+\$ \ °N 



Ajv-i, mi, m 2 + Ai, . . . , m N + Ajv_i 

(10) 



where 



N-l f) x ™ ( £n. 
1 TT " + 1 ^ PN 



with ^ = . . . , £n-i) and A n = X]fc=i^fc (defining A = 0). It is understood that the summation 
in (flOl) is (N — l)-fold, with the sum over the nth index £ n running from zero to L n — 1. This 
expression can be evaluated easily using mathematical software tools such as Mathematica. 

We may wish to consider the leading order expansion of P ,Nak at high SNR, which would 
yield expressions for the diversity and coding gains of the multihop link. To do this, we first note 
that the Meijer G-function used here is defined as a Mellin-Barnes integral where the integration 
path goes from — ioo to ioo such that it separates the poles of the integrand [[T6l §16.17]. It 
is straightfoward to show that this integral converges in our application and that the integrand, 
given by 

lNak{S) = Z{J) r(l - S) 

where 21(7) = {pNlth/ Y\ n &n) T - \ is well-behaved as s — > 00 in the right half s-plane. 
Moreover, I Nak (s) has poles at s = m n + A n _i + j for n — 0, . . . , N — 1 and j — 0, 1, . . ., and 
at s = when Atv_i = 0. Thus, we can employ the residue theorerro with the usual closing arc 

4 See, e.g., 1171 for an introduction to complex analysis and the residue theorem. 



DRAFT 



March 8, 2013 



in the right half plane to evaluate the G-functions in (flOl) at high SNR, which leads to a more 
accessible and intuitive asymptotic expression for the outage probability given bye 



p 



o,Nak 



rr (-)•>• ^2 ^2 ( c r _ i _ /)!/! 



"SCM*^)'^'^ <12) 

where m = min{m n }, ji denotes the multiplicity of m, and 

/ \ / \r ^Nak (s) 

u(s) = (s-m) — . 

z s 

The sets C%, . . . , £ M are disjoint sets of (N — l)-tuples {£i, . . . , £n-i} defined such that Jjv a fc is) 
has an rth order pole at s = m when {£i, . . . , € £ r . The details of the calculations that 

lead to this result (and a rigorous definition of the sets £ 1; . . . , £ M ) are given in the appendix. 

This expression can be evaluated easily for specific examples, but by retaining only the leading 
order term (i.e., r = /i), we arrive at the following general asymptotic equivalence: 

(log 7)^ 1 

Po,Nak ~ i>Nak — (13) 

where 

^•*-TC^"„ ^ }e£ ( _1) JTW\ (14) 

t«l)."!*JV-l JC^-M 

is the coding gain of the linlo. For the case where each hop fades independently of others and 
all hops experience nonidentically shaped fading (i.e., mi ^ ■ • ■ / m^), ;U = 1 and we have 

Po,Nak ~ (ygf E (-l)^' 1 ^^)^. (15) 

5 The notation g (7) — o(h (7)) signifies that lim^oo g (7) //i (7) = 0. 

6 This definition of the coding gain is slightly different to the one typically used in system analysis. In fact, the standard 
definition of the coding gain (see, e.g., [18]) cannot be applied here since it is only valid when Pout obeys a power law decay 
in 7. 
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In the other extreme where all hops experience identically shaped fading (i.e., m = rti\ 
m N ), jx = N and we can write 



-\N-1 



p (pNlth) m (log 7) .... 

r Q M a k ~ at • 1-tO) 

' (N - ly.mT (mf ilJZ 7 m 



The expressions given above demonstrate the well-known fact that mth order diversity is 
achieved in these cases. However, the analysis is useful in illustrating the rate (with respect 
to SNR growth) at which mth order diversity is attained. By taking the (usual) definition of 
diversity to be 

d = lim (17) 

7->oo — log 7 

we see from (fl"3l) that, for finite but large values of 7, 

d(j) = m -(„- i)™ + M\. ( i8) 



log7 \J0g7, 

Although this result points to slow convergence in terms of diversity (due to the log log 7/ log 7 
term), it does not illustrate the full picture since the coding gain must be taken into account. 
Indeed, through our analysis, we have presented an accurate expression for the coding gain of 
a fixed-gain AF multihop link, which to the best of the authors' knowledge has not yet been 
reported in the literature. The coding gain is, in general, a complicated expression, although 
it is straightforward to compute. However, we can draw some conclusions about the behavior 
of certain systems, such as those that experience identically shaped fading. For example, (fT6T ) 
points to the importance of having a well-designed destination receiver with a high-performance 
low-noise amplifier (LNA) when operating in the high SNR regime (pjv/7 = N 0j n must be as 
low as possible). 

Finally, it should be noted that the leading order results given here are not very accurate for 
low to mid-range 7 when /i > 1. This results from the fact that log 7 increases very slowly, 
which effectively means that all (log7) P 7~ m terms are of roughly the same order for finite 7 
and, thus, should be included in the approximation. In such a case, it is best to use the general 
expansion given by (PT2l) . 
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B. Weibull Fading 

If the channels adhere to a Weibull fading profile, the random variable X n has a Weibull density 
function. The corresponding PDF and Mellin transform are given in Table HI By substituting 
this transform into © and applying the definition of the Fox //-function lfl9ll along with the 
functional relations [1201 (2.8)] and [20, (2.11)], we can write the outage probability as 

(1,1) 



/ 



P ,Wei ~ I" ^2 ( _ 1) N 1 ^^hN+l 



a N 



(Ajv-i, 1) , ( 1 + Al , — ),..., 1 + — , — 



^n+l / Pn 



n=l x 



(19) 

where 

<pi= n 

n=l 

As with the Nakagami-m case, we can employ the residue theorem to obtain a simple asymp- 
totic expression for P 0) wei using elementary functions. Omitting the details (the methodology is 
the same as was described for Nakagami-m fading in the appendix), it is possible to derive the 
following asymptotic expression: 



X 

where again m = min {m n }, // denotes the multiplicity of m, and 

„ , ^ r r(A^ 1 - s )nl 1 r(i + ^ 
u{s) = {s- m ) fII -^ 

To leading order in 7, we have the asymptotic equivalence 



(log7) M 1 

Po,Wei ~ i>Wei — (21) 



where 

^ ,= w {<1 ' OT 

is the coding gain. Note that this analysis points to similar diversity and coding gain behavior 
as was discussed for Nakagami-m fading channels. 
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C. Rician and Hoyt Fading 

No closed form expression of P Q exists for the cases where all hops follow a Rician (or Hoyt) 
fading model. For the Rician case, we can employ the transform given in Table U along with 
Proposition [2] to obtain an expression for P ,mce m terms of a contour integral, which can be 
evaluated numerically using, for example, Mathematica. The residue theorem can be applied to 
this integral to derive an asymptotic expression in the form of 

p PNlt h Un(Kn + l) ^ . .. \Ajy_i * \P g^ 1 ^ (1) 

P °> mce TTf n ^ ^ ( } C ^(r-l -/)!/! 

x E (!) ( '°g — ft f - + - )'"" ^ + ° (?"') < 23 > 

p=0 



where 



and 



/V V 1LA / 



e (*) = - i) r r rn"l~) S) II r ( x + A "-i - s ) lFl + A «-i - s > ^ K «) • 

^ ' n=l 

To leading order in 7, we have 

N TS ,1/1 -\N-1 

p P7v7tfc tt + 1 (log 7) 

Po, ffice ~ -^yy 11 - g-^- • (24) 

Similarly, for the case where all hops follow a Hoyt fading model, we can follow the same 
procedure to derive the asymptotic expression 



^ ^ ^ (_) Xl ^{r-l-l)\l\ 



x t C) f tog PNlth u n }l» d) T ^ (25) 



where 



Or, \ !+ 2A iV-l N -1 /)A n / \ ^» / n_ \ l+2A n _i 
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and 

N 



f \ l 1 V r ( X n-i ~ s) n , \ j? 1 + A n _i - s 2 + A n _i - s f 1 - q 2 n \ 
71 {S) = {S ~ 1) T(l-s) J r ^ + n_1 - s ) ^ { 2 • 2 ' 1; VlTgl ) 

To leading order in 7, (1251) reduces to 

^ 1 , 2/1 -\N-1 

Po,^ ~ IF 3 TyI 11 - • (26) 

D. Semi-general Formula 

It is clear that similarities exist between the expressions for P Q given for the different fading 
distributions analyzed above. This observation leads to our second main result, which is in the 
form of a semi-general asymptotic formula for the outage probability of fixed-gain AF multihop 
links. 

Conjecture 3: Consider an iV-hop fixed-gain AF link, where the PDF of the channel power 
for each hop decays exponentially, and thus has a Mellin transform that is well-behaved at 
|Im (s)\ = 00. The outage probability of this link is given by 

li r-l n {r-\-l) , x 



v r- l-l)W 

r=l {i 1 ,...,i N _ 1 }ec r 1=0 v ' 



I 

X 

p=0 

as 7 —7- 00. In (|2~7T) . the constants A and are dependent upon the distributional parame- 



tQ(-l)'(>og^^ + o(7-"') (27) 



ters for the N channels in the link with B t being dependent upon the indices in the vector 
£ — (£1, . . . ,£n-i) as well; m defines the minimum shape parameter of the iV channel power 
distributions with \x being the multiplicity of m (e.g., m = 1 and fx = N if all hops are 
Rayleigh, Rician, or Hoyt); and C r (s) — (s — m) r I (s) (7 /A) s with / (s) being a simple form 
of the integrand in (J9) where all terms independent of s have been removed. 

This result is stated as a conjecture since we make no attempt to rigorously define the properties 
that the channel power PDFs must have to make it a theorem. Indeed, only a few channel 
models are typically employed in practice. These models include those discussed in Sections 
IIV-AI through IIV-CI We can easily apply the formula to those scenarios. 
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Additionally, this formula can be used to characterize the outage probability of inhomogeneous 
links. As a toy example, which is perhaps not often encountered in practice but serves as a useful 
illustration, consider a four-hop link where the per-hop channels are consecutively modelled as 
experiencing Nakagami-m, Weibull, Rician, and Hoyt fading. In this case, it is straightforward 
to show that the semi-general formula given by (l27t holds with the following definitions: 
PNlth (K 3 + 1) (1 + qjf 



A 



H EL ^ 

( 2g 4 \ 1+2A3 ^e + i lo, 



" r(mi)(^ 3 + l) Aa Vl + ?4 ' 11 



71=1 



C r (s) = (s - m) r r ^~ 1 f> T (mi -s)r(l+ 1 V (1 + A 2 - s) T (1 + A 3 - s) 



r (1 — s) V m 2 

lFl £ + A2 - S ' 15 ^ 2Fl ( 2 ' 2 5 1! \TTqJ 

m = min {mi,m,2, 1} . 
This example is discussed further in Section |V] 

E. Convergence of Asymptotic s 

It is natural to determine the conditions under which the leading order expressions for P Q given 
in the preceding section serve as good approximations to the actual outage probability. First, we 
study the question: how large does 7 have to be before "asymptotic" becomes "approximate"? For 
nonidentical Nakagami-m and Weibull fading scenarios, convergence occurs quickly since the 
leading order expression monotonically decreases with increasing 7. For quasi-identical fading, 
as well as Rician and Hoyt fading channels, we note that the asymptotic bound for P Q given 
above has the form P Q ~ b (log7) M ~ 1 7 _m where b is independent of 'y, fi > 2, and m > 1. 
For 7 > 1, it can be shown that this expression has a maximum at 7 — e^~ l ^ m , and is 
monotonically decreasing for 7 > e^~ l ^ m . Thus, it is necessary that 7 3> e^ _1 ^ m for the 
asymptotic expressions given above to be reasonable approximations for the outage probability. 
For hardened channels (i.e., when m ^> 1 for Nakagami-m and Weibull fading) or systems with 
only a few hops, this condition is satisfied easily in practice. 
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V. Simulation Results and Discussion 

In this section, we present numerical results obtained for homogeneous and inhomogeneous 
multihop links. For the homogeneous case, we simulated a number of scenarios in order to 
validate our analytical results. All four fading distributions mentioned above were considered, and 
different system parameters were chosen to illustrate the accuracy of our analytical framework 
as well as some interesting behavior of fixed-gain multihop links. For the inhomogeneous case, 
the four-hop example discussed in Section IIV-DI was studied. For all calculations, we assume 
only a first order asymptotic correction, i.e., L n = 2 for n = 1, . . . , N — 1. 

First, we present results for the case where each hop fades according to a Nakagami-m 
distribution. Fig. |2]through Fig. |4]illustrate the outage probability as a function of 7 for different 
three-hop links. Simulations are plotting along with analytical results obtained through the Meijer 
C7-function expression (fTOT t as well as the leading order asymptotic result dervied from the 
application of the residue theorem, which is given by ([T3l and (fl4l) . Furthermore, we compared 
our results to a lower bound on the probability of outage that was developed in [3|. This bound, 
which can be obtained through an application of a harmonic-geometric mean bound on the 
end-to-end SNR, has subsequently been applied and discussed in a number of works (e.g., HI). 

The accuracy of our analysis is clear from these examples. Moreover, although the asymp- 
totic expression converges to simulation results for large 7, we see that the Meijer G-function 
expression is a good approximation even at low and mid-range SNR. We also observe that the 
harmonic-geometric mean bound is loose in some cases, particularly when the inherent diversity 
in the channel decreases with each hop (i.e., m n decreases with increasing n). This behavior 
was also noted in flU. 

Fig. ID illustrates the effect that quasi-identical fading has on the convergence of the asymp- 
totic expression. In this example, the first and third hops yield the minimum shape param- 
eter (mi = m 3 = 1). This leads to a second order pole in the residue analysis, and thus 
Po.Nak ~ ?pNak (log 7) 7 -1 - The logarithmic term delays convergence as discussed in Section 
|IV-E[ and this is apparent in the figure in the plot of the leading order expansion. It is also 
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clear that the approximation can be made to be much more accurate for moderate SNR levels by 
refining the expansion through the inclusion of all correction terms of order (7 _1 ) (the curve 
labelled "Refined Asymptotic" in the figure]! 

Finally, we provide results for a five-hop system in Fig. |5] Here, we adopt the same shape 
factors that were used in [1 1 cf. Fig. 4 therein]; specifically, we have mi = m 2 = 5, m 3 = m 4 = 
2.5, and m 5 = 1.5. Again, there is excellent agreement between the simulation results and the 
Meijer G-function expression for P 0t Nak, an d the asymptotic result converges around 13 dB. In 
contrast, the harmonic-geometric mean lower bound diverges. 

In Fig. [6} we present results for homogeneous links where each hop experiences Weibull 
fading. Two, three, and four hop systems were considered where the shape parameters were 
chosen to be m\ = 1.5, m 2 = 2, m 3 = 2.5, and m 4 = 1. Again, the analytical results agree well 
with the simulations. Furthermore, we see how the addition of a fourth hop with less inherent 
diversity compared to the other hops affects performance. As previously mentioned, the hop with 
the minimum shape parameter dictates performance at high SNR, which is evident from the loss 
in diversity in the N = 4 hop link in the figure. 

Results for Rician and Hoyt channels are given in Fig. [7] and Fig. [H respectively. The 
degradation in performance with increasing numbers of hops is apparent for both cases. Moreover, 
the loss in finite SNR diversity can also be observed, which results from the log 7 terms in the 
numerator of the high SNR expansion for P a . The characterization of this behavior is beyond 
the scope of this paper, but is an interesting observation, nonetheless. 

Finally, we present results for an inhomogeneous link in Fig. |9j In this example, the fading 
channel corresponding to the first hop adheres to a Nakagami-m distribution with m\ = 2. 
The second hop channel follows a Weibull profile with m 2 = 1.5. The channels related to 
the third and fourth hops follow Rician and Hoyt distributions, respectively, with K 3 = 3 and 
#4 = 3/4. The scale parameter for the nth hop is given by 9 n = n/2 and we have defined 
p n = 1 — (n — 1) /10. In Fig. [91 the leading order asymptotic and the semi-general formula are 

7 The notation / (x) — Q (g (x)) implies 3 xo, k > such that / (a;) > k ■ g (x) for x > Xo. 
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plotted along with simulation results. Again, we see that although the graph points to convergence 
for the leading order expression in the asymptotic limit, the semi-general formula provides a 
much more accurate expression at low and mid-range SNR. Finally, it should be noted that this 
complex example cannot be studied with the theory that has been detailed in the literature to 
date, which exemplifies the versatility of our analytical framework. 



In this paper, we presented a novel, rigorous asymptotic analysis of the outage probability 
for fixed-gain AF multihop relay systems. Our analysis was general in nature, lending itself to 
application in a range of scenarios, including cases where the per hop fading processes adhere to 
completely different models and distributions. Specifically, we first provided a general asymptotic 
formula for the outage probability that is applicable to any system where the hops are statistically 
independent. We then provided analytical expressions for different fading distributions - namely 
Nakagami-m, Weibull, Rician, and Hoyt fading - and gave a brief discussion on the convergence 
of these formulae, which culminated in a semi-general closed-form formula for the outage 
probability at high SNR that can be applied to analyze homogeneous and inhomogeneous 
systems. Finally, we demonstrated through simulations that our theory is accurate, even at low 
to mid-range SNR in many cases of interest. 

Appendix A 
Properties of Mellin Transforms 

The Mellin transform of a real- valued function / (x) where x > is defined by 



VI. Conclusions 




and its inverse is given by 




(28) 
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for c > 0. A useful identity exists for the Mellin transform of derivatives of /, namely that Ifl4l 
6.1 (10)] 

M[fW;s]=(-l) t T ^^M\f;s-£\ (29) 

provided the first £ derivatives exist and are well behaved. We can also employ Mellin transforms 
to approximate the average of a function (with respect to the kernel x s_1 ) evaluated at a point 
close to, but to the right of x, i.e., f (x + e) for e > 0. This is done by expanding / near x, 
which yields 

/oo poo / L_1 J. \ L ~ l J- 

x^f (x + e)dx = j( x- 1 r£ (x) + O (e L ) J dx ~ £ \ M [/ W 5 s ] + ° ^) 

; (30) 



again, provided the derivatives existQ Note that this is an asymptotic expansion; convergence of 
the series is not guaranteed in general. 

If X is a random variable and fx [x) denotes its density function, then M [fx] s] = E [x* -1 ]. 
If fx (x) is defined for x > 0, then we can define the complementary cummulative distribution 
function (CCDF) 

POO 

Fx (x) = fx (t) dt. 

J X 

Now we have the following identity, which is calculated using integration by parts and, in 
particular, holds for probability distributions with exponentially decaying tailQ 

M[^x;s] = s- 1 M[f x ;s + l]. (31) 



Finally, suppose Z = Yln=i -^n where {X n } are statistically independent. Then it is easy to see 
that 



N 



M[f z ;s} = l[M[fx n ;s}. (32) 



n=l 



The notation ~ denotes asymptotic equivalence for e — > in this case. 

9 The identity is derived for the CCDF rather than the CDF in order to ensure the transform converges. 
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Appendix B 
Residue Calculations 

It is instructive to outline some of the residue calculations that were made to obtain the results 
given in this paper. Most calculations follow similar reasoning. Consequently, we only include 
calculations for Nakagami-m fading in this appendix. 

Consider the function I Nak (s) defined in (flTT) . Suppose there are T < iV unique shape 
parameters, where the tth parameter fn t has multiplicity fi t and Ylt=i = N. If Aat_i = 0, 



INak (s) has a simple pole at s = with residue res (Ijvafc, 0) = rin=i ^ ( m «)- Moreover, I 



Nak 



also has a pole at s = s q j = m q + A g _i + j for g G { 1 , . . . , AT} and j G {0, 1, . . .}, which is, in 

general, an rth order pole where 1 < r < N. The residue of I Nak at this pole is given by 

1 d r ~ x 
res (I Nak , s gJ ) = jT] {( s ~ s <?>i) r J ^« fc 



(r — 1)! s->s 9j 9s 

= -tt^wE ( r 7 'W 1 '" (s «> (log2)i **" (33) 

1 >' 1=0 V 7 

where the third equality follows from the Leibniz rule of differentiation of products (with z» (a) 
being the nth derivative of v evaluated at a) and u q j (s) = (s — s q j) r I^ak (s) -2 s - 

The residue theorem states that P D can be expressed as a series of residues of poles of I^ak 
IfTTl . Since 2 oc 7 _1 , if we wish to construct an approximation to P Q that is a function of the 
leading power of z (as 7 grows large), then we can ignore all residues for which j > 0. This 
leaves a finite summation of residues corresponding to different shape parameters. In fact, it is 
clear from (|33l ) that out of these residues, those at poles relating to the smallest shape parameter 
in the set {mi, . . . , m^} dominate for large 7. Thus, we only care about residues of the form 

res (I Nak , m q + A f _ x ) = -yy^yy E 7 ^ ^ ( ^ + Vl) (1 ° g Z ? 
where q E {1, . . . , N} such that m q < m n for all n, i.e., m q is the smallest shape parameter, 
which has multiplicity ft. In particular, we must determine conditions under which A§_i = 
since this yields the leading order in z. 
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Denote the ordered indices of the corresponding minimum shape parameters by qi, . . .,qp,. 
Now, since {Ai, . . . , Ajv-i} are cummulative sums of the indices {£1, . . . , £n-i}, it follows that 
In ok (s) will have a fith order pole at m q when Aq A _i = 0, a (ft — l)th order pole at m q when 
Ag A = but A^„_i > 0, and so on. We wish to enumerate these instances through the indices 
{£i, . . . ,£at_i}. To this end, it is possible to construct (x disjoint sets of (N — l)-tuples, which 
we denote C\,..., Cp_, such that C r consists of all sets {£i, . . . , £n-i} that satisfy the conditions 
Aq r -i = and Aq r+1 _i > 0, where the second condition is only necessary (and valid) for 
r < p, — 1. It follows that the order of the pole at s = rriq is r if an only if {£±, . . . , £n-i} € C r . 
This allows us to partition the expression of P a into terms related to the order of the poles. 
Finally, we can substitute the leading order residues, summing over the appropriate sets of £ 
indices, to write the following expression for P Q 

l In \ nj r =i {e 1 ,...,e N _ 1 }ec r 1=0 K ' 

Sustituting for z and applying the binomial theorem yields the result given in (fT2l) . 
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R1 V- -H R2 
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Fig. 1. Multihop system diagram with one source transmitting to a destination via N — l relay nodes; no direct source-destination 
link exists. 




Fig. 2. Po.Nak vs. 7 for a fixed-gain AF multihop systems with Nakagami-ra fading channels (N = 3; m n = n, K n — 2 and 
On = Pn = 1 for all n). 
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Fig. 3. P ,Nak vs. 7 for a fixed-gain AF multihop systems with Nakagami-m fading channels (N — 3; m„ = N — n + 1, 
K n — 2, and 8„ — p n — 1 for all n). 



TABLE I 

Mellin transforms for various fading distributions. 



Distribution 


Channel Power Density Function fx(x), x > 


Mellin Transform M {fx; s + 1} 


Nakagami-m 

Weibull 

Rician 

Hoyt 


6»-™T(m)' 1 a: m - 1 e- 3;/9 

m Q-m x m-l e -{x/e) m 

e- 1 (k + 1) e -{K+o-HK+v*) Io (Jm^mA 
i±d e ^r?f-' Io (k=sL x \ 

2q9 c 10 \ I^e x ) 


8 s r (s + m) /r (m) 

6> s r (s/m+ 1) 

e ~ K (l?TT) Sr ( s + 1 ) xM' + Mi*) 
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Fig. 4. P ,Nak vs. 7 for a fixed-gain AF multihop systems with Nakagami-m fading channels (N — 3, mi = 7713 = 1, 
m 2 =2, pi = 1, p 2 = 1/3, p 3 = 5/3; and 0„ = (N - n+ 1) /2 and AT„ = 2 for all n). 
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Fig. 5. Po.Nak vs. 7 for a fixed-gain AF multihop systems with Nakagami-m fading channels (N = 5, mi = m2 = 5, 
?Ti3 = m4 = 2.5, ?Ti5 = 1.5; and K n — 2 and 8 n = p n = 1 for all n). 
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Fig. 6. P ,Wei vs - 7 for a fixed-gain AF multihop system with Weibull distributed channels (JV = 2, 3, 4, mi = 1.5, m,2 = 2, 
7713 = 2.5, = 1; and A" n = 2 and # n = p n = 1 for all n). 
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7 (dB) 



Fig. 7. Po,mce vs. 7 for a fixed-gain AF multihop system with Rician distributed channels (N = 2, 3, 4, ATi = 1, K2 = 3, 
-K3 = 5, K4 — 0; and K„ = 2 and # n = p n = 1 for all n). 
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Fig. 8. P ,Ho V t vs. 7 for a fixed-gain AF multihop system with Hoyt distributed channels (N — 2, 3, 4, qi = 3/4, q-z — 1/2, 
qa = 1/3, qi = 1/4; and if n = 2 and (9 n = p n = 1 for all n). 



DRAFT 



March 8, 2013 



29 



k. . . yt ! ! ! 


X Simulated 




........... , 0< ............................... 

: %* 

'v X 


Leading Order Asymptotic 

- i - Semi-general Formula 




x 

N X 






■V v 




r '-^ „ % ~ - 




vx 

/ X V 




N x 

x 


X 

X' 


" x 

. X 

X 


V xX 


X 

X 


X 

nX 

X X 

V :: >X - 


% ■■■■\ - 




X x X * 

- x x - 


.. .. V V.: \ ^ 











5 10 15 20 25 30 35 

7 (dB) 



Fig. 9. P vs. 7 for a four-hop inhomogeneous fixed-gain AF multihop system. The channels corresponding to the four hops 
follow a Nakagami-m, Weibull, Rician, and Hoyt distribution, respectively with the following parameters: mi = 2, m,2 = 1.5, 
K 3 = 3, g 4 = 3/4, n = n/2, and p n = 1 - (n - 1) /10. 
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